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Abstract 

Using the effective potential approach for composite operators, we have for- 
mulated a general method of calculation of the truly nonperturbative Yang- 
Mills vacuum energy density in the covariant gauge QCD ground state quan- 
tum models. It is defined as an integration of the truly nonperturbative part 
of the full gluon propagator over the deep infrared region (soft momentum 
region). A nontrivial minimization procedure makes it possible to determine 
the value of the soft cutoff in terms of the corresponding nonperturbative 
scale parameter, which is inevitably present in any nonperturbative model 
for the full gluon propagator. We have shown for specific models of the full 
gluon propagator explicitly that the use of the infrared enhanced and finite 
gluon propagators lead to the vacuum energy density which is finite, always 
negative and it has no imaginary part (stable vacuum) , while the infrared van- 
ishing propagators lead to unstable vacuum and therefore they are physically 
unacceptable. 

PACS numbers: ll.15.Tk, 12.38.Lg 
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I. INTRODUCTION 



The nonperturbative QCD vacuum is a very complicated medium and its dynamical and 
topological complexity [1-3] means that its structure can be organized at various levels (clas- 
sical, quantum). It can contain many different components and ingredients which contribute 
to the truly nonperturbative vacuum energy density (VED), one of the main characteristics 
of the QCD ground state. Many models of the QCD vacuum involve some extra classical 
color field configurations such as randomly oriented domains of constant color magnetic 
fields, background gauge fields, averaged over spin and color, stochastic colored background 
fields, etc. (see Refs. [1,4,5] and references therein). The most elaborated classical models 
are random and interacting instanton liquid models (RILM and IILM, respectively) of the 
QCD vacuum [6]. These models are based on the existence of the topologically nontrivial 
instanton-type fluctuations of gluon fields, which are nonperturbative, weak coupling solu- 
tions to the classical equations of motion in Euclidean space (see Ref. [6] and references 
therein). 

Here we are going to discuss the quantum part of VED which is determined by the 
effective potential approach for composite operators [7-9]. It allows us to investigate the 
nonperturbative QCD vacuum, in particular Yang-Mills (YM) one, by substituting some 
physically well-justified Ansatz for the full gluon propagator since the exact solutions are not 
known. In the absence of external sources the effective potential is nothing but VED which 
is given in the form of the loop expansion where the number of the vacuum loops (consisting 
in general of the confining quarks and nonperturbative gluons properly regularized with the 
help of ghosts ) is equal to the power of the Plank constant, h. 

The full dynamical information of any quantum gauge field theory such as QCD is 
contained in the corresponding quantum equations of motion, the so-called Schwinger- 
Dyson (SD) equations for lower (propagators) and higher (vertices and kernels) Green's 
functions. It is a highly nonlinear, strongly coupled system of four-dimensional inte- 
gral equations for the above-mentioned quantities. The kernels of these integral equa- 
tions are determined by the infinite series of the corresponding skeleton diagrams [10-12]. 
It is a general feature of nonlinear systems that the number of exact solutions (if any) 
can not be fixed a priori. Thus formally it may have several exact solutions. These equa- 
tions should be also complemented by the corresponding Slavnov- Taylor (ST) identities 
[10-12] which in general relate the above mentioned lower and higher Green's functions to 
each other. These identities are consequences of the exact gauge invariance and therefore 
"are exact constraints on any solution to QCD" [10]. Precisely this system of equations 
can serve as an adequate and effective tool for the nonperturbative approach to QCD [13,14]. 

Among the above-mentioned Green's functions, the two-point Green's function describing 
the full gluon propagator 

iD^{q) = {T,Md(-q\i) + iL^(q)} ^, (1-1) 

has a central place [10-15]. Here £ is a gauge fixing parameter (£ = 0, Landau gauge) 
and T^ u {q) = g^ u — q^q v /q 2 = g^ u — L^ u (q). Evidently, its free perturbative (tree level) 
counterpart is obtained by simply setting the full gluon form factor d(—q 2 , £) = 1 in Eq. 
(1.1). In particular, the solutions of the above-mentioned SD equation for the full gluon 



2 



propagator (1.1), are supposed to reflect the complexity of the quantum structure of the QCD 
ground state. As emphasized above, it is a highly nonlinear system of four-dimensional 
integrals containing many different, unknown in general, propagators, vertices and kernels 
[10-12]. Because of truncation schemes this system becomes the equation for the full gluon 
propagator only, but it remains nonlinear, nevertheless. Different truncations could lead to 
qualitatively different solutions, and the number of these solutions may be increased only. 
Moreover, to clearly distinguish between the exact or approximate solutions (if any), we do 
not know even the complete set of boundary conditions to attempt to uniquely fix solution 
of the truncated equation. We certainly know the boundary condition in the ultraviolet 
(UV) limit because of asymptotic freedom and certainly we do not know the corresponding 
boundary condition in the infrared (IR) because precisely of confinement (at this stage it 
is not even clear whether the two boundary conditions (in the UV and in the IR (if it 
can be established)) will be sufficient to completely fix the theory or not). Because of the 
above-discussed highly complicated mathematical structure of the SD equation for the full 
gluon propagator, there is no hope for exact solution(s). However, in any case the solutions 
of this equation can be distinguished by their behaviour in the deep IR limit (the UV 
limit is uniquely determined by asymptotic freedom), describing thus many different types 
of quantum excitations and fluctuations of gluon field configurations in the QCD vacuum. 
Evidently, not all of them reflect the real structure of the QCD vacuum. 

The deep IR asymptotics of the full gluon propagator can be generally classified into the 
three different types: 1) the IR enhanced (IRE) or IR singular (IRS), 2) the IR finite (IRF) 
and 3) the IR vanishing (IRV) ones (for references see the corresponding sections below). Let 
us emphasize that any deviation in the behaviour of the full gluon propagator in the deep 
IR domain from the free perturbative one automatically assumes its dependence on a scale 
parameter (at least one) in general different from QCD asymptotic scale parameter, Aq C d- 
It can be considered as responsible for the nonperturbative dynamics (in the IR region) in 
the QCD vacuum models under consideration. If QCD itself is a confining theory, then such 
characteristic scale is very likely to exist. In what follows, let us denote it, say An p. This 
is very similar to asymptotic freedom which requires the above-mentioned asymptotic scale 
parameter associated with nontrivial perturbative dynamics in the UV region (scale viola- 
tion). However, for calculation of the truly nonperturbative VED we need not exactly the 
deep IR asymptotics of the full gluon propagator, but rather its truly nonperturbative part, 
which vanishes when the above-mentioned nonperturbative scale parameter goes formally to 
zero, i.e., when only the perturbative phase survives. So we define the truly nonperturbative 
part of the full gluon form factor in Eq. (1.1) as follows: 

d NP (-q\ A NP ) = d(-q 2 , A NP ) - d(-q 2 , A NP = 0), (1.2) 

which, on one hand, uniquely determines the truly nonperturbative part of the full gluon 
propagator. On the other hand, the definition (1.2) explains the difference between the truly 
nonperturbative part d NP (—q 2 ) and the full gluon propagator d(—q 2 ) which is nonperturba- 
tive itself. Let us note in advance, that in the realistic models for the full gluon propagator, 
the limit A NP — > is usually equivalent to the limit — q 2 — > oo. In some cases, the model 
gluon propagator does not depend explicitly on the nonperturbative scale parameter (the 
dependence is hidden) then its behaviour at infinity should be subtracted. In the realistic 
models of the full gluon propagator its truly nonperturbative part usually coincides with its 
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deep IR asymptotics underlying thus the strong intrinsic influence of the IR properties of 
the theory on its nonperturbative dynamics. 

It is well known, however, that VED in general is badly divergent in quantum field theory, 
in particularly QCD [16]. Thus the main problem is how to extract the truly nonperturbative 
VED which is relevant for the QCD vacuum quantum model under consideration. It should 
be finite, negative and it should have no imaginary part (stable vacuum). Why is it so 
important to calculate it from first principles? i.e., on the basis of some realistic Ansatz 
for the full gluon propagator only. As was emphasized above, this quantity is important 
in its own right as being nothing else but the bag constant (the so-called bag pressure) 
apart from the sign, by definition [16]. Through the trace anomaly relation [17] it helps in 
the correct estimation of such an important phenomenological nonperturbative parameter 
as the gluon condensate introduced in the QCD sum rules approach to resonance physics 
[18]. Furthermore, YM VED assists in the resolution of the U(l) problem [19] via the 
Witten-Veneziano (WV) formula for the mass of the meson [20]. The problem is that 
the topological susceptibility [19-23] needed for this purpose is determined by a two point 
correlation function from which perturbative contribution is already subtracted, by definition 
[20,23-25]. The same is valid for the above-mentioned bag constant which is much more 
general quantity than the string tension because it is relevant for light quarks as well. Thus 
to calculate correctly the truly nonperturbative VED means to understand correctly the 
structure of the QCD vacuum in different models. 

We have already formulated a general method of calculation of the truly nonperturbative 
YM VED in the axial gauge QCD in Ref. [26], where the Abelian Higgs model [27] of the 
dual QCD [28] ground state was investigated. Moreover, we have calculated the truly non- 
perturbative VED (using particular method) in the covariant gauge QCD quantum vacuum 
model as well [29,30]. The main purpose of this paper (section II) is to formulate precisely a 
general method of calculation of the truly nonperturbative quantum part of YM VED in the 
covariant gauge QCD. In sections III, IV and V it is illustrated by by considering different 
covariant gauge QCD quantum models of its ground state by choosing three different types 
of the deep IR asymptotics of the full gluon propagator, IRE, IRF and IRV, respectively. 
The conclusions are presented in section VI. 

II. THE TRULY NONPERTURBATIVE VACUUM ENERGY DENSITY 

In this section we formulate a general method of numerical calculation of the quantum 
part of the truly nonperturbative YM VED in the covariant gauge QCD. Let us start from 
the gluon part of VED which to-leading order (log-loop level ~ fr)[] is given by the effective 
potential for composite operators [7] as follows: 

V{D) = \\ ^ Tr[HD ° lD) ~ {D ° lD) + 1} ' (2 ' 1} 



1 Next-to- leading and higher contributions (two and more vacuum loops) are numerically sup- 
pressed by one order of magnitude in powers of h at least and are left for consideration elsewhere. 
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where D(q) is the full gluon propagator (1.1) and D (q) is its free perturbative (tree level) 
counterpart. Here and below the traces over space-time and color group indices are under- 
stood. The effective potential is normalized as V(D ) = 0, i.e., the free perturbative vacuum 
is normalized to zero. In order to evaluate the effective potential (2.1) we use the well-known 
expression 



TrlniD^D) = 8 x In det^ 1 D) = 8 x 4 In 



(2.2) 



It becomes zero (in accordance with the above mentioned normalization condition) when the 
full gluon form factor is replaced by its free perturbative counterpart. This composition does 
not depend explicitly on a gauge choice. Going over to four-dimensional (n = 4) Euclidean 
space in Eq. (2.1), on account of (2.2), and evaluating some numerical factors, one obtains 
(e 9 = V(D)) 



~hS dq 



V 



ln[l + 3d(q 2 )] - -d(q 2 ) + a 



(2.3) 



where constant a = (3/4) — 2 In 2 = —0.6363 and the integration from zero to infinity 
is assumed. Substituting the definition (1.2) into the Eq. (2.3) and doing some trivial 
rearrangement, one obtains 



7T Z 



dq 2 q 2 



ln[l + 3d^(q 2 , A NP )} - -d NP (q 2 ,A NP ) 



7T Z 



IPT, 



where we introduce the following notation 



IpT 



2 2 

Q 



ln[l + 



3d(q 2 ,A NP = 0) 
l + 3d NP (q\A NP ) 



3 

-d(q 2 ,A NP = 0) + a 



(2.4) 



(2.5) 



It contains the contribution which is mainly determined by the perturbative part of the full 
gluon propagator, d(q 2 ,A]yp = 0). The constant a also should be included since it comes 
from the normalization of the free perturbative vacuum to zero. If we separate the deep IR 
region from the perturbative one (which consists of the intermediate (IM) and UV regions 
since the IM region remains terra incognita in QCD) by introducing the so-called soft cutoff 
explicitly then we get 



~ / ° da 

7T Z JO 



2 2 

Q 



\n[l + 3d NP (q 2 ,A NP )} - ^d NP (q 2 ,A NP ) 



— (I PT + I PT ), (2.6) 



where evidently 

„ roo 

Ipt = dq 

J 1o 



V 



ln[l + 3d NP (q 2 , A NP )} - U NP (q 2 , A NP ) 



(2.7) 



Thus the first integral represents contribution to YM VED which is determined by the truly 
nonperturbative piece of the full gluon propagator integrated over the deep IR region. In 
other words, just this term is the truly nonperturbative contribution to YM VED. This 
means that the two remaining terms in Eq. (2.6) should be subtracted by introducing 
corresponding counter terms into the effective potential. Thus in general the integral (2.5) 
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V{G) = -ij ^-TriHG^G) - (G l G) + 1}, (2.8) 



determining the contribution from the perturbative part of the full gluon propagator and the 
integral (2.7) determining the contribution from the perturbative region (IM plus UV) are of 
no importance for our present consideration. The above-mentioned necessary subtractions 
can be done in more sophisticated way by means of ghost degrees of freedom (see below). 
The effective potential at the log-loop level for the ghost degrees of freedom is 

d n p 

(27T)"" 

where G{p) is the full ghost propagator and Go(p) is its free perturbative (tree level) coun- 
terpart. The effective potential V(G) is normalized as V(G ) = 0. Evaluating formally the 
ghost term e g h = V(G) in Eq. (2.8), we obtain e g h = n~ 2 I gh . The integral I g h depends on the 
ghost propagator, which remains arbitrary (unknown) within our approach. In principle, we 
have to sum up all contributions to obtain the total VED (the confining quark part of the 
vacuum energy density is not considered here). However, upon the substitution of definition 
(1.2) into the integral over the whole momentum range from zero to infinity (2.3), some 
terms appear there which may have unphysical singularities below the scale Aqcd (integral 
(2.5)). Thus the initial VED (2.3) is formal one, it suffers from unphysical singularities 
briefly mentioned above and it is badly divergent as well. In order to get physically mean- 
ingful expression, one have to subtract two integrals (2.5) and (2.7) from Eq. (2.3). We 
have done this subtraction with the help of ghost term by imposing the following condition 
A = Ipt + Ipt — Igh = 0. The nonperturbative gluon contribution to VED is determined by 
subtracting unwanted terms by means of the ghost contribution, i.e., defining e g + e g h = £ym 
at A = 0. Thus the truly nonperturbative YM VED becomes 



7T 



1 f q , 2 2 

"2 / d<i q 
T 2 JO 



d NF {q\ A NP ) - ln[l + M NF {q\ A NP )} 



(2.9) 



In many cases this subtraction is sufficient to obtain the expression for the truly nonpertur- 
bative YM VED. However, in some other cases the truly nonperturbative part of the full 
gluon propagator which enters Eq. (2.9) continues to suffer from unphysical singularities 
below the scale Aqcd (see discussion at the end of section V). As it was noticed, some 
additional terms should be included in our subtraction scheme in this case indicating that 
the chosen Ansatz for the full gluon propagator itself was not realistic. 

A few general remarks are in order. In QCD nothing should explicitly depend on ghosts. 
By contributing into the closed loops only, the main purpose of their introduction is to cancel 
unphysical degrees of freedom of gauge bosons (maintaining thus the unitarity of S- matrix), 
for example to exclude the longitudinal components, the above mentioned unphysical sin- 
gularities below the QCD scale, etc. This is the main reason why they are to be considered 
together with gluons always. In nonperturbative QCD in general and in our approach in 
particular the ghost propagator (or equivalently the ghost self-energy) still remains unknown 
(in this sense arbitrary) since the exact ghost-gluon vertex (which enters the corresponding 
SD equation) is not exactly known (in Refs. [31,32] some very specific truncation scheme 
is used in order to derive particular expression for this vertex). We know however, that 
the ghost propagator contribution to VED, regular or singular, should be combined with 
gluon contribution in order to cancel exactly the above-mentioned unphysical singularities 
of the gauge bosons which are inevitably present in any Ansatz for the full gluon propaga- 
tor. In other words, if one knows the ghost propagator exactly then the above-mentioned 
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cancellation should proceed automatically (as usual in perturbative calculus if, of course, 
all calculations were correct). But if it is not known exactly (as usual in nonperturbative 
calculus) then one has to impose condition of cancellation as it was done in our case, A = 0. 
Obviously, the above-mentioned condition of cancellation was imposed in the most general 
form. Instead of the introduction of some counter terms into the initial effective potential 
to cancel the most dangerous UV divergences presented in the integral (2.5), we have used 
the ghost term for this purpose as well. Thus our subtraction scheme is in agreement with 
general physical interpretation of ghosts to cancel all unphysical degrees of freedom of the 
gauge bosons [10,33]. 

The expression (2.9) is our definition of the truly nonperturbative YM VED as integrated 
out the truly nonperturbative part of the full gluon propagator over the deep IR region (soft 
momentum region, < q 2 < q$). The soft cutoff (as a function of the nonperturbative 
scale) can be determined by the corresponding minimization procedure (see below). 



A. 



From this point it is convenient to factorize scale dependence of the truly nonperturbative 
YM VED (2.9). As was already emphasized above, d NP (q 2 ) always contains at least one 
scale parameter (A^p) responsible for the nonperturbative dynamics in the model under 
consideration. It is considered as free one within our general method, i.e., "running" (when it 
formally goes to zero then the perturbative phase only survives in the model). Its numerical 
value (if any) will be used at final stage only to evaluate numerically the corresponding 
truly nonperturbative YM VED (if any). We can introduce dimensionless variables and 
parameters by using completely extra scale (which is always fixed in comparison with Atvp), 
for example flavorless QCD asymptotic scale parameter A Y m as follows: 

Zo = ^-' b= w- (210) 

•M.i/ IV YM YM 

Here z is a corresponding dimensionless soft cutoff while the parameter b has a very clear 
physical meaning. It measures the ratio between nonperturbative dynamics, symbolized by 
A 2 NP and nontrivial perturbative dynamics (violation of scale, asymptotic freedom) symbol- 
ized by Ay M . When it is zero only perturbative phase remains in the model. In this case, 
the gluon form factor obviously becomes a function of z and b, i.e., d NP (q 2 ) = d NP (z, b) and 
the truly nonperturbative VED (2.9) is {eyu = £ym(^o,&)) 

flg(z , b) = —^—e Y M(zo, b), (2.11) 

A VM 

where the gluon effective potential at a fixed scale, A Y m, [26,29,34] is introduced 



1 f z o 

Q g = Q g (z , b) = — dz z 



^d NP (z, b) - m[l + 3d NP (z, b)} 



(2-12) 



This expression precisely allows us to investigate the dynamical structure of the YM vacuum. 
It is free of scale dependence since it has been already factorized in Eq. (2.11). It depends 
only on z and b and a minimization procedure with respect to b, dfl g (z , b)/db = 0, (usually 
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after integrated out in Eq. (2.12)) can provide a self-consistent relation between zq and 6, 
that is we get go as a function of A^p. Let us note in advance that the final numerical 
results will depend on A NP only as it should be for the nonperturbative part of YM VED 
(see sections III and IV below). Obviously, the minimization with respect to z leads to 
trivial zero. In principle, through the relation A YM = q^ZQ 2 , it is possible to fix soft cutoff 
qo itself, but this is not the case indeed since then zq can not be varied. 



B. 



On the other hand, the scale dependence can be factorized as follows: 

„2 



Z = 



A 2 ' 
IV NP 



A 2 ' 



(2.13) 



i.e., 6=1. For simplicity (but not loosing generality) we use the same notations for the 
dimensionless set of variables and parameters as in Eq. (2.10). In this case, the gluon form 
factor obviously becomes a function of z only, d NP (q 2 ) = d NP (z) and the truly nonpertur- 
bative YM VED (2.9) becomes 



_ 4 -2 

2% Z 



zo 



dz z 



d NP (z) - ln[l + 3d 



NP, 



(2.14) 



Evidently, to fix the scale is possible in the two different ways. In principle, we can fix A^p 
itself, i.e., introducing 



tt g (z ) 



1 



A 4 

NP 



1 [ z ° 

n z Jo 



dz z 



d NP (z) - ln[l + 3d NP {z)} 



(2.15) 



However, the minimization procedure again leads to the trivial zero, which shows that this 
scale can not be fixed. 

In contrast to the previous case, let us fix the soft cutoff itself, i.e., setting [26,29,30] 



- 1 1 2 f Z " 

qo TV jo 



dz z 



^d NP (z) - ln[l + 3d NP (z)] 



C2-16) 



In this case the perturbative phase is recovered in the z — > oo (A NP — > 0) limit. Now the 
minimization procedure with respect to z is nontrivial. Indeed, dCl g (z ) / dz Q = 0, yields the 
following "stationary" condition 



dz z 



-d NP (z) - ln[l + 3d NP (z)] 



l -z 2 



^d NP (z )-Hl + 3d NP (z )} 



[2.17) 



the solutions of which (if any) allow one to find qo as a function of A^p. On account of this 
"stationary" condition, the effective potential (2.16) itself becomes simpler for numerical 
calculations, namely 



2tt 2 



d NP {4) - ln[l + 3d NP {4 t )} 



(2.18) 



where Zq is a solution (if any) of the "stationary" condition (2.17) and corresponds to the 
minimum(s) (if any) of the effective potential (2.16). In the next sections, we illustrate how 
this method works by considering some quantum models of the covariant gauge QCD ground 
state explicitly. 
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III. THE IRE GLUON PROPAGATOR. ZME QUANTUM MODEL 



Today there are no doubts left that the dynamical mechanisms of the important non- 
perturbative quantum phenomena such as quark confinement and dynamical (or equivalently 
spontaneous) chiral symmetry breaking (DCSB) are closely related to the complicated topo- 
logically nontrivial structure of the QCD vacuum [1-4,10]. On the other hand, it also becomes 
clear that the nonperturbative IR dynamical singularities, closely related to the nontrivial 
vacuum structure, play an important role in the large distance behaviour of QCD [35,36]. For 
this reason, any correct nonperturbative model of quark confinement and DCSB necessarily 
turns out to be a model of the true QCD vacuum and the other way around. 

Our model of the true QCD ground state is based on the existence and importance 
of such kind of the nonperturbative, quantum excitations of the gluon field configurations 
(due to self-interaction of massless gluons only, i.e., without explicit involving some extra 
degrees of freedom) which can be effectively correctly described by the g~ 4 behaviour of the 
full gluon propagator in the deep IR domain (at small q 2 ) [29,30]. These excitations are 
topologically nontrivial also since they lead to the nontrivial YM VED (see below). Thus 
our main definition (1.2) becomes 

d NP (-q 2 , A NP ) = d(-q 2 , A NP ) - d(-q 2 , A NP = 0) = (3.1) 

(-<r) 

In the above-mentioned papers [29,30] the nonperturbative scale was denoted as p,, i.e., 
p, = Ajvp. In this way we obtain the generally accepted form of the deep IR singular 
asymptotics for the full gluon propagator (for some references see below) 

D,M)-{q 2 )~\ ? 2 -0, (3.2) 

which may be refered equivalently to as the strong coupling regime [10]. It describes the zero 
momentum modes enhancement (ZMME) dynamical effect in QCD at large distances. We 
prefer to use simply ZME (zero modes enhancement) since we work always in momentum 
space. This is our primary dynamical assumption in this section. The main problem due to 
this strong singularity is its correct treatment by the dimensional regularization method [37] 
within the distribution theory [38], which was one of highlights of our previous publications 
[29,30] (see also Ref. [39]). There exist many arguments in favor of this behaviour: 

a) . Such singular behaviour of the full gluon propagator in the IR domain leads to the 
area law for static quarks (indicative of confinement) within the Wilson loop approach [40]. 

b) . The cluster property of the Wightman functions in QCD fails and this allows such 
singular behaviour like (3.2) for the full gluon propagator in the deep IR domain [41]. 

c) . After the pioneering papers of Mandelstam in the covariant (Landau) gauge [42] 
and Baker, Ball and Zachariasen in the axial gauge [43], the consistency of the singular 
asymptotics (3.2) with direct solution of the SD equation for the full gluon propagator in 
the IR domain was repeatedly confirmed (see for example Refs. [13,14,44,45] and references 
therein). 

d) . Moreover, let us underline that without this component in the decomposition of the 
full gluon propagator in continuum theory it is impossible to "see" linearly rising potential 
between heavy quarks by lattice QCD simulations [46] not involving some extra (besides 
gluons and quarks) degrees of freedom. This should be considered as a strong lattice evidence 
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(though not direct) of the existence and importance of q~ 4 -tjpe excitations of gluon field 
configurations in the QCD vacuum. There exists also direct lattice evidence that the zero 
modes are enhanced in the full gluon propagator indeed [47]. 

e) . Within the distribution theory [38] the structure of the nonperturbative IR singu- 
larities in four-dimensional Euclidean QCD is the same as in two-dimensional QCD, which 
confines quarks at least in the large N c limit [48]. In this connection, let us note that 
q~ 4 IR singularity is the simplest nonperturbative power singularity in four- dimensional 
QCD as well as q~ 2 IR singularity is the simplest nonperturbative power singularity in 
two-dimensional QCD. The QCD vacuum is much more complicated medium than its two- 
dimensional model, nevertheless, the above-mentioned analogy is promising even in the case 
of the nonperturbative dynamics of light quarks. 

f) . Some classical models of the QCD vacuum also invoke q~ 4 behaviour of the gluon 
fields in the IR domain. For example, it appears in the QCD vacuum as a condensation of the 
color-magnetic monopoles (QCD vacuum is a chromomagnetic superconductor) proposed by 
Nambu, Mandelstam and 't Hooft and developed by Nair and Rosenzweig (see Ref. [49] and 
references therein. For recent developments in this model see Di Giacomo's contribution in 
Ref. [1]) as well as in the classical mechanism of the confining medium [50] and in effective 
theory for the QCD vacuum proposed in Ref. [51]. 

g) . It is also required to derive the heavy quark potential within the recently proposed 
exact renormalization group flow equations approach [52]. 

h) . It has been shown in our papers that the singular behaviour (3.2) is related directly 
to light quarks confinement and DCSB [29,30]. Moreover, a very good agreement has been 
obtained with the phenomenological values of the topological susceptibility, the mass of the 
r]' meson and the gluon condensate [21,22]. 

Thus we consider our main Ansatz (3.1), (3.2) as physically well-motivated. Let us 
emphasize that d NP (—q 2 ,^) = A 2 NP /(—q 2 ) is the truly nonperturbative part of the full gluon 
propagator since it vanishes in the perturbative limit (A^ P — ► 0, when the perturbative 
phase survives only) and simultaneously it correctly reproduces the deep IR asymptotics of 
the full gluon propagator, i.e., d NP (—q 2 ) coincides with d IR (—q 2 ). 



A. 

The truly nonperturbative YM VED is given now by Eq. (2.9) with d NP (q 2 ) = A 2 NP /q 2 
in Euclidean space. Let us first introduce the A-type set of dimensionless variables (2.10). 
Then d NP (q 2 ) becomes d NP (z, b) = b/z. Performing almost trivial integration in the effective 
potential at a fixed scale (2.12), one obtains 



Q g (z ,b) 



1 



2tt 2 



9 ^ |1 + I)-t 



b — z"l In 



1 + 



3b' 



(3.3) 



It is easy to show that as a function of b, the effective potential (3.3) linearly approaches 
zero from below and it diverges also linearly at infinity while as a function of z itself it 
approaches zero from above and also diverges as ~ — Zq at infinity. Thus as a function of b it 
has a local minimum (relating b to z ) at which the truly nonperturbative YM VED will be 
always finite and negative. The minimization procedure with respect to b, dfl g (z ; b)/db = 0, 
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yields the following "stationary" condition, v = 41n(l + (z//3)), where v = zo/b. Its solution 
is u mm = 2.2. Using this "stationary" condition, the effective potential (3.3) can be written 
down as follows: 



2 , .rain 



b 2 V 



2tt 2 



In 1 + 



so the truly nonperturbative YM VED (2.11) becomes 



= -0.12736 2 , 



(3.4) 



e YM = -0.1273 x A 



(3.5) 



where the relation A^ P = b 2 Ay M has been already used. Determined in this way, it is 
always finite (since characteristic scale of our model A^p is finite, evidently it can not 
be arbitrary large), automatically negative (as it should be for the truly nonperturbative 
energy) and it has no imaginary part (stable vacuum). Obviously the characteristic scale 
of our model A^p can not be determined within the YM theory alone. Its numerical value 
should be taken from good physical observable in full QCD by implementing the physically 
well-motivated scale-setting scheme. Precisely this has been done in our papers [29,30] 
where the nonperturbative VED was numerically evaluated from first principles. Moreover, 
in recent publications [21,22] it is shown that our numerical results are of the necessary order 
of magnitude in order to nicely saturate the large mass of rj' meson in the chiral limit as well 
as the phenomenological value of the topological susceptibility. Thus the existence of the 
nontrivial VED in ZME quantum model, which agrees well with QCD topology, is one more 
serious argument in its favor. It is worthwhile to present numerical value for the soft cutoff 
in terms of A NP , namely q = 1.48324Ajvp- It follows from the solution of the "stationary" 
condition, of course. 



B. 

It is instructive to calculate the truly nonperturbative YM VED by choosing the B-type 



set of dimensionless variables (2.13). Then d (q 



A 2 NP /q 2 becomes d 



NPi 



l/z. 



Performing almost trivial integration in the effective potential at a fixed scale (2.16) in this 
case, one obtains 



2tt 2 



9 In 1 + 



z%ln 



1 + 



(3.6) 



It is easy to show now that as a function of z , the effective potential (3.7) diverges as ~ Zq 1 
at small z and converges as ~ — Zq 1 at infinity (perturbative limit), see Fig. 1. Thus as 
a function of z it has a local minimum at z = 41n(l + (z /3)), the so-called "stationary" 
condition in this case. Its solution again is z™ m = 2.2. At "stationary" state the effective 
potential (3.6) can be written down as follows: 



O / min\ 



2tt 2 



r)- 1 - In ( 1 + -t- 

\ Z 



-0.0263, 



(3.7) 



so the truly nonperturbative YM VED (2.16) becomes 
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€YM — 



-0.0263g^ = 



-0.1273 x At 



(3.8) 



where the relation q$ = (z™ m ) 2 A^ p has been already used. Thus we have explicitly demon- 
strated that truly nonperturbative YM VED does not indeed depend on how one intro- 
duces dimensionless variables into the effective potential, i.e., €ym = A^r P r2 s (z/ mm , b) = 
qQQ g (z™ m ) = — 0.1273A^ P . In some cases, the B-type calculation is preferable. For ex- 
ample, to calculate confining quark contribution into the total VED is much easier using 
precisely this set of the dimensionless variables (see our papers [29,30] and next section as 
well). 
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FIG. 1. Effective potential (3.7) as a function of zq. 



IV. THE IRF GLUON PROPAGATOR 

Let us consider now a possible IRF behaviour of the full gluon propagator (in the Landau 
gauge) in the deep IR domain, which was suggested by recent lattice calculations in Ref. 
[53]. The main definition (1.2) in this case becomes 

d»?(- q \ M) = d(- q \ M) - d(- q \ M = 0) = ^l^-fi - (4.1) 

Here M is the mass scale parameter responsible for the nonperturbative dynamics in this 
model, i. e., M = A^p in our notation. When the parameter M formally goes to zero, 
the perturbative phase only remains in this model. Again as in previous case, the truly 
nonperturbative part vanishes in the perturbative limit (M — > 0) and it reproduces the 
IR asymptotics of the full gluon propagator correctly as well. The best estimates for the 
parameters M and A are M = (1020 ± 100 ± 25) MeV and A = (9.8 + 0.1 - 0.9) . As it was 
emphasized above, the numerical value of the parameter M will be used only at final stage 
in order to estimate numerically the truly nonperturbative YM VED in this model. The 
exponent in general is a = 2 + 5, where 5 > and small, while Z « 1.2 is the renormalization 
constant. 
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In this case, it is convenient to choose the B-type set of variables and parameters (2.13). 
Then d NP (q 2 ) in Euclidean space becomes 

dNP ( z ) = 7i — ~~~vn~> ( 4 - 2 ) 

v ' (1 + z) 1+a v ' 

where the parameter a\ = ZA = 11.76 is fixed. Substituting this into the effective potential 
(2.16), one obtains 

1 1 

Cl g (z ; a{) = -^t Y M = x z 2 {li(z ; ai) - I 2 {z ; ai)}, (4.3) 
where integrals are given as follows: 

SdiZ 



h(z ;a 1 )=J dzzln (l + °\* 1+a ), 



(l + z) ] 

3a Z ° z 
I 2 (z ; ai ) = ~yJ dzz (1+ Z z) i +a - ( 4 - 4 ) 
o v ; 

The asymptotic behaviour of the effective potential (4.3) depends on the asymptotic proper- 
ties of the integral I\{zq] a\) since the integral ^{zq] a\) in Eq. (4.4) can be taken explicitly, 
namely (in what follows in this section, a — 2) 

Hz ; oi) = ^(ln(l + z ) + 2[(1 + zo)- 1 - 1] - + ^o)- 2 - 1])- (4.5) 

From these expressions it is almost obvious that asymptotics of the effective potential (4.3) at 
z — > 0, oo to-leading order can be easily evaluated analytically. The effective potential (4.3) 
as a function of the soft cutoff zo has two local minimums (see below). The corresponding 
"stationary" condition can be evaluated as follows: 

[/,(„;«,) - /,(»,«,)] = iz5{jn(l + ^) - jJ^j}. (4.6) 

Using this "stationary" condition, the effective potential (4.3) at "stationary" state becomes 

ft,^)--^!*^)- a^U, (4.7) 



where Zq* is a solution(s) to the "stationary" condition (4.6). The two solutions of the 
"stationary" condition (4.6) corresponding to the two local minimums are Zq = 0.19 and 
Zq = 2.37 with almost equal numerical values for the corresponding effective potentials at 
"stationary" states, namely fi 9 (0.19;ai) = —0.0309 and f2 g (2.37;ai) = —0.0310, respec- 
tively. However, the numerical values of the nonperturbative YM VED (4.3) are drastically 
different, 

e yM (0.19) = -0.0309go(0.19) = -0.00123 x M 4 (4.8) 

and 

e yM (2.37) = -0.0310g 4 (2.37) = -0.174 x M 4 , (4.9) 

where the relation = (z^ m ) 2 M 4 and the corresponding values of z™ n (= Zq) were applied. 
How to distinguish between the two solutions for the truly nonperturbative YM VED (4.8) 
and (4.9)? This question is discussed in the following. 
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A. Discussion 



In the first case, on account of the numerical value of the nonperturbative scale M pa 
1 GeV, Eq. (4.8) numerically becomes 

e yM (0.19) = -0.00123 GeV 4 . (4.10) 

It is the same order of magnitude as VED due to instantons [22]. Thus summing up this 
and instantons with ZME values, one obtains a fair agreement with chiral QCD topology 
[20]. Also the soft cutoff in this case is qo ~ 0.463M pa 463 MeV. This is quite reasonable 
value for the deep IR region (in continuum theory) where the smooth-type behaviour of the 
full gluon propagator effectively takes place. 

In the second case, on account of the numerical value of the nonperturbative scale M pa 
1 GeV, Eq. (4.9) numerically becomes 

e yM (2.37) = -0.174 GeV 4 . (4.11) 

In Refs. [21,22] an analytical formalism has been developed which allows one to calculate 
the topological susceptibility as a function of the truly nonperturbative YM VED. The 
corresponding expression is 




(4.12) 



where parameter £ has two different values, namely ^,nsvz = 2/11 and ^ HZ = 4/33 (see 
Ref. [22]). Evaluating (4.12) numerically, on account of (4.11), one obtains, xf SVZ — 
(550.8 MeV) 4 and X ? Z = (259.6 MeV) 4 , while its phenomenological value is, xt™ = 
(180.36 MeV) 4 . Thus, Eq. (4.11) substantially overestimates the phenomenological value of 
the topological susceptibility (in both modes) and consequently the mass of rj' meson in the 
chiral limit, indeed. The soft cutoff in this case is go ~ 1.54M ~ 1.54 GeV. It is also hard 
to imagine that the deep IR region (in continuum theory) can be effectively extended up 
to pa 1.54 GeV especially for the smooth-type behaviour of the full gluon propagator there. 
The continuum limit of the scale parameter M is not known, so its realistic numerical value 
still remains to be well-established, and so does the selection from solutions, Eqs. (4.8) and 
(4.9). Let us note, that in accordance with the general scheme of our method we distinguish 
the nonperturbative scale of this model from the perturbative one but for simplicity we 
retain the same notation. Evidently, one will obtain the same numerical results for the truly 
nonperturbative YM VED by choosing the set of variables of A type. 



V. THE IRV GLUON PROPAGATOR 



The IRV full gluon propagator is represented by the so-called Zwanziger-Stingle (ZS) 
formula [54,55] 
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in the whole range, where /i 4 is again the mass scale parameter responsible for the non- 
perturbative dynamics in this model, i.e., fi = A^p, in our notation. When it is zero 
then the ZS gluon propagator (5.1) becomes free perturbative one, indeed. Though the full 
gluon propagator (5.1) is nonperturbative itself, however its truly nonperturbative part is 
determined by the subtraction (1.2), i.e, 

d NP {q 2 ) = d(q\ ^) - d(q 2 , / = 0) = - 2 ( . (5.2) 

Since this expression is rather simple, it will be instructive to perform calculations in both 
schemes, A (2.10) and B (2.13). So let us start from A scheme. 

A. 

Within the A-type set of variables (2.10), d NP (q 2 ) from Eq. (5.2) becomes d NP (z,b) = 

— (b 2 /b 2 + z 2 ) (Euclidean space). After the integration over four-dimensional Euclidean space 
in Eq. (2.12), one obtains 

n g (z , b) = ^{~8b 2 \n2b 2 + 8b 2 \n(2b 2 - z 2 ) + (b 2 + Az 2 ) \n(b 2 + z 2 ) 

-4z 2 Hz 2 -2b 2 ) -b 2 \nb 2 }. (5.3) 

From this expression it follows obviously that the effective potential (5.3) at any finite 
relation between the soft cutoff z and parameter b will always contain the imaginary 
part, which is a direct manifestation of the vacuum instability [56] in this model. Its 
asymptotics at b — > 0, oo to-leading order can be easily evaluated analytically. Omit- 
ting all intermediate calculations, one finally obtains, Q g (z ,b) — (9/87r 2 )6 2 ln6 2 and 
Q g (zo,b) ~f>^oo — (l/8' 7r2 )[3 + 4 ln(— 2)]^ 2 , confirming the vacuum instability. Let us also 
consider the corresponding formal "stationary" condition, dfl g (zo,b)/db = 0, which yields 

3t 2 + (1 + t 2 ) ln(l + tl) + 8(1 + t 2 ) ln(l - |) = 0, (5.4) 

where tjj = (z 2 /b 2 ). It has only trivial solution to = Zo = 0. 

Thus the vacuum of this model is unstable, indeed, so it has no relation to quark con- 
finement and DCSB. Our conclusion is in full agreement with conclusion driven in Ref. [57]. 
The particular-type expressions for the dressed-quark-gluon vertex free from ghost contri- 
butions were used in their investigation. Our result, however, is a general one since it does 
not require the particular choice of the dressed-quark-gluon vertex. 

B. 

Within the B-type set of variables (2.13), d NP (q 2 ) from Eq. (5.2) becomes d NP (z) = 

— (1/1 + z 2 ) (Euclidean space). After almost trivial integration over four-dimensional Eu- 
clidean space in Eq. (2.16), one obtains 
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^(*o) = ^o~ 2 {- 8 In 2 + 8 ln(2 - z 2 ) + (1 + 4z 2 ) ln(l + z 2 ) - 4z 2 ln(z 2 - 2)}. (5.5) 

From this expression it obviously follows that the effective potential at any finite value of 
the soft cutoff Zq will always contain the imaginary part, which is a direct manifestation 
of the vacuum instability [56] as it was indicated above. Its asymptotics at z Q — > 0, oo to- 
leading order can be easily evaluated analytically. Omitting all intermediate calculations, 
one finally obtains, 0, g (z ) ~ zo -+o — (l/87r 2 )[3 + 41n(— 2)] and Cl g (z ) ~ 20 _ >O o (9/87r 2 )^ 2 ln^, 
so the vacuum of this model is unstable, indeed. In order to confirm this, let us consider 
the corresponding formal "stationary" condition which is 

3z 2 + (1 + zl) ln(l + zl) + 8(1 + z 2 ) ln(l - |) = 0. (5.6) 

It has only trivial solution z = 0. 

In Ref. [57] it was proposed the modification of the ZS propagator (5.1) which took into 
the consideration the renormgroup improvements to-leading order for the running coupling 
constant in the UV region, namely 

= x . . (5.7) 



V A QCD J 



Here const obviously depends on the first coefficient of the /^-function and an unphysical 
parameter r is introduced in order to regulate unphysical singularity - Landau pole - at 
q 2 = Aq CD (Euclidean space). The truly nonperturbative part now is 

i'V) = d(q\ /) - d(,V = 0) = - , £ , 4 x , C ° nS \ x . (5.8) 



V ^QCD J 



However, it is possible to show that YM VED continues to contain imaginary part in this 
case as well. It is worth noting, that in derivation of the corresponding expression for YM 
VED (2.9) all terms depending in general on some unphysical parameters (in this case r) 
should be additionally subtracted by means of ghosts (as it was mentioned above in Section 
II just after Eq. (2.9)). Concluding, let us note that neither (5.2) nor (5.8) coincides with 
deep IR asymptotics of the corresponding full gluon propagators (5.1) and (5.7). 



VI. CONCLUSIONS 



In summary, we have formulated a general method how to numerically calculate the 
quantum part of the truly nonperturbative YM VED (the bag constant, apart from the 
sign, by definition) in the covariant gauge QCD quantum models of its ground state using 
the effective potential approach for composite operators. It is defined as integrated out 
the truly nonperturbative part of the full gluon propagator over the deep IR region (soft 
momentum region), Eq. (2.9). The nontrivial minimization procedure makes it possible to 
determine the value of the soft cutoff as a function of the corresponding nonperturbative scale 
parameter which is inevitably present in any nonperturbative full gluon propagator model. 
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If the chosen Ansatz for the full gluon propagator is realistic one, then our general method 
gives the truly nonperturbative YM VED which is always finite, automatically negative and 
it has no imaginary part (stable vacuum) (sections III and IV). Its numerical value does not, 
of course, depend on how one introduces the scale dependence by choosing different scale 
parameters as it was described above in subsections A and B of section II, i.e., both set of 
variables lead to the same numerical value of the truly nonperturbative YM VED. 

From comparison of Eqs. (2.3) and (2.9), a prescription can be derived how one can 
obtain the relevant expression for the truly nonperturbative YM VED. For this purpose the 
full gluon propagator in Eq. (2.3) should be replaced by its truly nonperturbative part in 
accordance with Eq. (1.2). The constant a should be omitted (it was already explained 
why) and the soft cutoff q$ on the upper limit should be introduced. Now it looks like the 
UV cutoff. Nevertheless let us underline once more that it separates the deep IR region from 
the perturbative one, which includes the IM region as well. It has a clear physical meaning 
as determining the range where the deep IR asymptotics of the full gluon propagator is 
valid. By definition it can not be arbitrary large as the UV cutoff is. As far as one chooses 
Ansatz for the full gluon propagator, the separation "NP versus PT" is exact because of 
the definition (1.2). The separation "soft versus hard" momenta is also exact because of the 
above-mentioned minimization procedure. Thus the proposed determination of the truly 
nonperturbative YM VED is uniquely defined. The nontrivial minimization procedure can 
be done only by two ways. First, to minimize the effective potential at a fixed scale (2.11), 
(2.12) with respect to the physically meaningful parameter. When it is zero, the perturbative 
phase only survives in all models of the QCD ground state. Equivalently, we can minimize 
the auxiliary effective potential (2.16) as a function of the soft cutoff zq itself. When it 
goes to infinity then again the perturbative phase survives only. On the other hand, both 
effective potentials (2.12) and (2.16) should go to zero in the perturbative limit since the 
perturbative contributions have been already subtracted from the very beginning (see section 
2). As was underlined above, both methods lead to the same numerical value for the truly 
nonperturbative YM VED. 

We have shown explicitly that the IRE gluon propagator (3.2) as well as IRF (4.1) 
correspond to nontrivial VED which is always finite, negative and it has no imaginary part 
(stable vacuum). In this way they reflect some physical types of excitations of gluon field 
configurations in the QCD vacuum. At the same time, the IRV gluon propagators (5.1) and 
(5.2) lead to unstable vacuum and therefore are physically impossible. However, these results 
are by no means general. For example, to come to the same conclusion for the IRV gluon 
propagator obtained and investigated in Refs. [31,32] it is necessary to proceed along with 
lines of our method. Thus the proposed method is precisely general one and each particular 
model for the full gluon propagator should be separately analyzed within its framework. 
However, it seems to us that the unstable vacuum is a fundamental defect of all vacuum 
models based on the IRV-type behaviour of the full gluon propagator. It is worthwhile also 
noting that, in contrast to IRE gluon propagator, the smooth behaviour of the full gluon 
propagator in the IR domain is hard to relate to quark confinement and DCSB. 

Thus our method can serve as a test of any different QCD vacuum models (quantum 
or classical) since it provides an exact criterion for the separation "stable versus unstable 
vacuum" . Vacuum stability in classical models is important as well. For example, we have 
already shown [26] that the vacuum of the Abelian Higgs model without string contributions 
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is unstable against quantum corrections. 

There is no general method of calculation of the confining quark contribution to the total 
VED. In quantum theory it heavily depends on the particular solutions of the corresponding 
quark SD equation, on account of the chosen Ansatz for the full gluon propagator. If it 
is correctly calculated then it is of opposite sign to the nonperturbative gluon part and it 
is one order of magnitude less (see, for example our papers [21,22,29,30]). Our method is 
not a solution for the fundamental badly divergent problem of VED in QCD. Moreover, it 
is even not necessary to deal with this problem. What is necessary indeed, is to be able 
to extract the finite part of the truly nonperturbative VED in a self-consistent way. Just 
this is provided by our method which thus can be applied to any nontrivial QCD vacuum 
quantum/classical models. 

In conclusion, let us make some remarks. In some cases together with the nonpertur- 
bative scale some other parameter (s) should be considered as "running" in accordance with 
the general scheme of our method. For example, such situation will arise in the IRF model 
gluon propagator suggested by lattice calculations in Ref. [58] (see also Ref. [59]). In this 
case the general procedure of calculation of the truly nonperturbative YM VED (if any) 
remains, of course, unchanged. However, due to some technical details (for example, the 
corresponding "stationary" condition (2.17) will be more complicated) this case requires a 
separate consideration. A brief recent reviews on both continuum and lattice gluon prop- 
agators can be found in Refs. [15,53]. An attempt of VED calculation by introduction of 
rather controversial gluon mass was made in a recent paper [60]. 
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